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Abstract 


By analyzing large-aspect-ratio spiral-defect-chaos (SDC) convection images, 
we show that the Karhunen-Loeve decomposition (KLD) scales extensively for 
subsystem-sizes larger than Ad (d is the fluid depth), which strongly suggests 
that SDC is extensively chaotic. From this extensive scaling, the intensive 
length ^kld is computed and found to have a different dependence on the 
Rayleigh number than the two-point correlation length Local computa¬ 
tions of ^kld reveal a substantial spatial nonuniformity of SDC that extends 
over radii 18d < r < 45d in a F = 109 aspect-ratio cell. 
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A significant theoretical challenge is to hnd ways to characterize the nonperiodic time- 
dependent patterns often observed in large sustained nonequilibrium systems In a re¬ 
cent paper |^, a new length scale for characterizing spatiotemporal chaos (STC), the KLD 
length ^kld, was proposed based on the extensive scaling of the Karhunen-Loeve decomposi¬ 
tion (KLD) ip. This length scale was shown to be computable from moderate amounts of 
space-time data and to contain information similar to the fractal-dimension density, which 
has not yet been computed from experimental data [^]. For certain idealized mathemat¬ 
ical models, the length .^kld was shown to have a different parametric dependence than 
the commonly computed two-point correlation length ^2 and so provides a different way to 
characterize spatiotemporal chaos 0]. Further, this length .^kld could be calculated from 
data localized to a region of space and so offered a way to analyze spatial inhomogeneities. 
However, an application of the KLD length to experimental data had not been made prior 
to this work. 

In this paper, we present the hrst application of the KLD length to experimental STC 
data by analyzing the recently discovered ^ spiral-defect chaos (SDC) in Rayleigh-Benard 
convection. We analyze the SDC state of large aspect-ratio (F = r/d = 29 and 109 where r 
and d are the radius and thickness of the cell) cells and hnd that the KLD dimension Dkld 
of the data scales extensively with subsystem volume when the diameter of the subvolume 
is larger than about Ad where d is the depth of the huid. This extensivity strongly suggests 
that the SDC state is extensively chaotic which provides the hrst such evidence for an 
experimental system. From the extensive scaling of the dimension Dkld with subvolume, 
we calculate the length ^kld which rehects the density of linearly independent modes needed 
to approximate the spatiotemporal data 0. For the reduced Rayleigh numbers e explored 
in this study, ^kld exhibits a minimum with increasing e while the two-point correlation 
length ^2 monotonically decreases. We speculate that this diherent behavior in the KLD 
length arises from changes in the structure of spirals and straight rolls with increasing e. 
Finally, we also show that the length .^kld provides a way to measure a dynamical inhomo¬ 
geneity in SDC data within the interior of the cell. A local KLD analysis demonstrates that 
SDC is not dynamically uniform (Most numerical simulations of SDC [0 have used periodic 
lateral boundary conditions, for which the dynamics should be statistically homogeneous by 
translational invariance). In addition, this nonuniformity extends to within a radius r = ISd 
of the center of an aspect-ratio F = 109 cell. These results suggest that the length ^kld can 
provide useful insights about spatiotemporal chaos. 

We tested the utility of the length ^kld for analyzing SDC experimental data |^,|^ by 
collecting thousands of shadowgraph images 0 of the SDC state from a Rayleigh-Benard 
experiment using compressed CO 2 gases with Prandtl number a ~ 1 and in cylindrical cells 
of aspect ratio F = 29 and 109. SDC appeared for e = AT/ATc — 1 above 0.56 and 0.23 
respectively for the F = 29 and 109 cells, where ATc is the critical temperature difference 
for the onset of convection. Five sets of SDC data were taken in the two cells as summarized 
in Table |. We chose sampling times At that were signihcantly longer than the correlation 
time r ~ lOt^ (where t^ is the vertical thermal diffusion time) to maximize the amount of 
uncorrelated data 0 . Good spatial resolution of the convection patterns required at least 5 
pixels per convection roll for the F = 109 cell. Data was taken from the central 43% of the 
F = 109 convection cell, as indicated in Fig. |l]. 

To illustrate the KLD analysis used to study the experimental data, let u(ti,yij) denote 
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the light intensity of an experimental image at position at time tj. The analysis proceeds 
by constructing a T x S' space-time matrix of data, 


Aij = < u(ti,yij) > 


( 1 ) 


where < u{ti, Xj) > denotes the time-average (average over index i) of the set of experimental 
images u{ti, xj), T is the number of observation times ti, and S is the number of observation 
sites Xj. The KLD dimension -Dkld |13 matrix A^j then measures the number of 

linear eigenmodes needed to approximate some fraction 0 < / < 1 of the variance of the 
experimental data and can be computed from the eigenvalues of the matrix A^A 


We compute T)KLD(xj) for concentric subsystems of volume V (square or circular geometry) 
that are centered at a particular point xj in space. Zi)KLD(xj) depends on the point xj 
and so provides a measure of dynamical inhomogeneity. If Zi)KLD(xj) increases linearly with 
subsystem volume V and with a slope 5, then the length .^kld is dehned to be where 
d is the dimensionality of the data (d = 2 for SDC). 

In applying these ideas to SDC data, we computed Dkld for a hxed fraction / = 0.7 
and for larger and larger square subimages in the center of the convection cell of size S x S, 
where 2d < S < 25d. As shown in Fig. ^ the dimension Dkld scales approximately linearly 
with subsystem data over a relatively large range of subsystem sizes Ad < S < 13d provided 
that a sufficiently long time series was used. This extensive scaling, together with the 
arguments in Ref 0 relating .^kld to the dimension correlation length ^s, strongly suggests 
that the SDC state is extensively chaotic. We note that the extensive linear scaling of 
the dimension Dkld with subsystem area is best for smaller subsystems, which we believe 
is a consequence of the fact that smaller subsystems have a faster time scale to become 
statistically stationary. 

We next compared ^kld (computed in the center of the cell for / = 0.7) with the 
two-point correlation length ^2 as the reduced Rayleigh number e was varied. The two- 
point correlation length ^2 was calculated from the inverse of the width of the peak in the 
Fourier spectrum of the spatial data which was pre-multiplied with a hanning window of 
diameter equal to the lateral dimension of the images. As shown in Table ||, the parametric 
dependences of the length ^kld and the two-point correlation length ^2 are different |T2 . 
The fact that .^kld attains a minimum between e = 0.79 and 0.88, and then increases 
with increasing e (corresponding to a decrease in complexity) is somewhat counterintuitive 
since one might have expected .^kld to decrease with e as the system is forced further 
away from equilibrium (more modes are needed to approximate the space-time data). A 
possible explanation for these opposing trends may be found by examining the different 
spatial structures in Figs. HA and P For Fig. |A (e = 0.52), the area fraction of local 
straight-roll regions is larger than that for spiral-roll regions, whereas for Fig. (e = 0.93) 
the relationship between straight- and spiral-roll regions is reversed. We speculate that the 
data between e = 0.79 and 0.88 consist of nearly equal fractions of straight- and spiral-roll 
regions and thus require more KLD eigenmodes per unit volume DKLo/fsub and so .^kld = 
(DKLD/Vsub)~^A is smaller. This speculation about the relative complexity of straight- 
roll versus spiral-roll regions is supported by examining the KLD spatial eigenmodes for 
the e = 0.52 and e = 0.93 data. These modes have the same qualitative symmetries up 
to eigenmode 21 (/=0.45 for e = 0.52), beyond which the local straight-roll regions in the 
e = 0.52 data entered into the next KLD spatial modes and broadened the KLD eigenvalue 
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spectruml^. In the e = 0.93 data, beyond eigenmode 21 the eigenmodes consisted of lattices 
of convection rolls. 

As the KLD analysis of a snbsystem is a local procedure, one can quantify differences in 
the SDC in different regions of the convection cell and so test the assumption that SDC is 
dynamically homogeneous. For the F = 109 cell, we investigated a radial dependence of .^kld 
by hrst establishing that Hkld scaled extensively for an annulus of radius Vo < r < Vo + d 
and angular sector that varied from A6 = 0 to A6 = radians (extensivity was with 
respect to AQ). We could not estimate nonuniformity for < Sd because the subsystems 
were too small, or for > 45d as this exceeded the area imaged by the CCD camera. 
Fig. ^ shows how the length .^kld increases by 10% from r = 8d to r = 45d. The variation 
in the length ^kld with radial distance r demonstrates that SDC cannot be considered 
homogeneous (Fig. even in large aspect ratio cells. Based on this calculation we consider 
the cell approximately uniform for r < 18d. The time-average of the SDC data does not 
resemble the nonuniformity in the KLD length in the interior of the cell, but the more easily 
computed variance pattern does resemble that of ^kld- The 10% nonuniformity of ^kld in 
the large cell may be due to both dynamical and experimental reasons. The dynamical one 
would be due to the pinning of the pattern by the side-wall which has been shown to cause 
time-averaged patterns near the sidewall for instance in rotating convection experiments . 
(We also found nonunformity of ^kld near the sidewall in the F = 29 cell.) The experimental 
source could come from the nonlinearity of shadowgraph method. For thin cells and high 
epsilon values, the nonlinearity is strong and can make the shadowgraph method sensitive 
to small optical nonuniformities. Unfortunately, the relative strength of these two sources 
of nonuniformity cannot be determined. 

In summary, we have carried out the first analysis of experimental data using the ex¬ 
tensive scaling of KLD in subsystems for the spiral-defect chaos state. KLD analysis is 
straightforward to apply to small subsystems of different geometry as it does not impose an 
approximate periodicity of the space-time data as is the case for Fourier analysis. By verify¬ 
ing that the dimension Dkld scaled linearly with subsystem size, we provide strong evidence 
that the experimental system is extensively chaotic [|l|. The utilization of a subsystem was 
essential for our study for two reasons. First, local analysis of subsystems allowed the char¬ 
acterization of inhomogeneous dynamics in our experimental convection cell. Further, by 
exploiting the reduced data requirements of subsystems, we could estimate the parametric 
behavior of the dimension density (the average number of degrees of freedom per unit area) 
in the center of the convection cell. For the SDC data analyzed, the length ^kld seems to 
quantify differences in the fraction of straight-roll and spiral-defect regions (also observed in 
the KLD spatial eigenmodes), and so provides information beyond that available from ^ 2 - 
Finally the nonuniformity of the KLD length was shown to extend over 80% of the radius 
in the F = 109 convection cell. Future analysis will hopefully provide further insights about 
the transition to SDC and about the chiral-symmetry breaking recently observed in SDC 
for rotating convection experiments |jl4 . 
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FIGURES 


FIG. 1. Time-dependent shadowgraph patterns for the central 43% of the T = 109 convection 
cell, with light and dark regions corresponding to cool and warm fluid respectively. (A) and (B) 
are snap-shots of the SDC pattern for reduced Rayleigh number e = 0.52 and e = 0.93 respectively. 


FIG. 2. Scaling of the KLD dimension Ukld with subsystem area N'^ for data in the central 43% 
of the r = 109 cell with reduced Rayleigh number e = 0.52 and variance fraction / = 0.70. The 
area of a subsystem is measured in units of d where d is the the thickness of the cell. The labels 
indicate the number of images used in the calculations and lines connecting points were drawn to 
guide the eye. 

FIG. 3. Radial inhomogeneities of the length .^kld in annular subsystems of fixed radii and 
increasing azimuthal angle for the F = 109 convection cell (e = 0.79). The radial distance r 
is measured in units of d (the thickness of the convection cell) from the center of the F = 109 
convection cell. 
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TABLES 

TABLE 1. Parameters for data used to cal¬ 
culate ^kld- L is the aspect ratio, e is the re¬ 
duced Rayleigh number, T is the number of im¬ 
ages, and At/tv is the sampling rate in units of 
the vertical thermal diffusion time tv 


r 

e 

T 

At/tv 

tv (s) 

29 

1.80 

2048 

33 

8.28 

109 

0.52 

3100 

1050 

0.86 

109 

0.67 

2500 

1050 

0.86 

109 

0.79 

2500 

1050 

0.86 

109 

0.93 

2570 

728 

0.82 


TABLE II. Lengths ^kld and ^2 (normalized 
to the depth of the fluid d) as functions of aspect 
ratio r and reduced Rayleigh number e. The 
length ^kld was calculated for fraction f = 0.7. 


r 

e 

^KLo/d 

6/d 

29 

1.80 

1.47 

2.28 

109 

0.52 

0.84 

4.38 

109 

0.67 

0.81 

3.85 

109 

0.79 

0.75 

3.60 

109 

0.88 

0.74 

3.39 

109 

0.93 

0.90 

3.22 
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